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Abstract 

The relationship between observed variability time and emission region ge- 
ometry is explored for the case of emission by relativistic jets. The approximate 
formula for the jet-frame size of the emission region, R' = DcAt b s is shown to 
lead to large systematic errors when used together with observed luminosity and 
assumed or estimated Doppler factor D to estimate the jet-frame photon energy 
density. These results have implications for AGN models in which low-energy 
photons are targets for interaction of high energy particles and photons, e.g. 
synchrotron-self Compton models and hadronic blazar models, as well as mod- 
els of intra-day variable sources in which the photon energy density imposes a 
brightness temperature limit through Compton scattering. 

The actual relationship between emission region geometry and observed vari- 
ability is discussed for a variety of geometries including cylinders, spheroids, 
bent, helical and conical jet structures, and intrinsic variability models including 
shock excitation. The effects of time delays due to finite particle acceleration 
and radiation time scales are also discussed. 
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1 Introduction 

In the standard picture of active galactic nuclei (AGN), accretion onto a super-massive 
black hole is via an accretion disk, and a significant fraction of the accretion power 
(possibly supplemented by tapping into the rotational energy of the black hole) pro- 
duces twin opposing relativistic jets moving outward along the disk axis, with typical 
Lorentz factors T ~ 2-10 as inferred from very long baseline interferometry (VLBI) 
observations. The objects observed in high energy 7 rays are "blazars", AGN in which 
one of the jets is closely aligned toward the observer. It is natural that in 7-rays we 
should see preferentially those AGN with aligned jets because the emission from the 
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jet is Doppler boosted in energy and relativistically beamed along the jet direction (for 
a discussion of relativistic effects see Urry and Padovani 1995). The 7 ray emission 
from blazars is variable (as it is also at optical, UV and X-ray energies). Relativistic 
effects also cause the observed variability time to be shorter than the time scale over 
which the emission changes in the jet frame. 

The spectral energy distribution (SED) of blazars shows two broad peaks, the low 
energy peak extending from the infrared to the UV or X-ray region of the spectrum, 
and the high energy peak starting in the X-ray or 7 ray range. The usual interpretation 
is that relativistic electrons produce the low energy part by synchrotron emission, and 
that the same electrons produce the high energy part by Compton scattering the low 
energy part and/or external photons to higher energies. The 3rd EGRET catalog of 
high-energy 7-ray sources (Hartman et al. 1999) contains around 70 high confidence 
identifications of AGN, and all appear to be blazars (Montigny et al. 1995, Mukherjee 
et al. 1997). Clearly, the 7-ray emission is associated with AGN jets. 

Four BL Lac objects have been detected in the TeV energy range: Mrk 421 (Punch 
et al. 1992), Mrk 501 (Quinn et al. 1996), IE S2344+514 (Catanese et al. 1998) and 
PKS 2155-304 (Chadwick et al. 1999). Recently, the spectrum of Mrk 501 has been 
measured up to 24 TeV by the HEGRA telescopes (Konopelko et al. 1999). Several 
of the EGRET AGN show 7-ray variability with time scales of ~ 1 day (Kniffen et 
al. 1993) at GeV energies. The TeV 7-ray emission of two BL Lacs shows very rapid 
variability. For Mrk 421, variability on a time scale as short as ~ 15 minutes has been 
reported (Gaidos et al. 1996). In the case of Mrk 501, variability on a time scale of a few 
hours was observed during the 1997 high level of activity, and there is evidence of a 23 
day periodicity (Protheroe et al. 1998, Hayashida et al. 1998) interpreted in terms of a 
binary black hole model for the central engine by Rieger and Mannheim (2000). These 
variability timescales place important constraints on the models. For example, the 
synchrotron self-Compton (SSC) model appears to be just consistent with recent multi- 
wavelength observations of Mrk 421 and Mrk 501 during flaring activity (Bednarek and 
Protheroe 1997, 1999). However, the allowed range of physical parameters (Doppler 
factor and magnetic field) is rather small, and this mechanism may well be excluded 
by future observations. For a recent review of TeV 7-ray astronomy see Kifune (2002). 

Rapid variability in intra-day variability (IDV) sources is a long-standing problem 
as it implies apparent brightness temperatures in the radio regime which may exceed 
10 17 K or relativistic beaming with extremely high Doppler factors, coherent radia- 
tion mechanisms, or special geometric effects (Wagner and Witzel 1995). The very 
rapid flaring observed at TeV and X-ray energies during flaring activity in blazars also 
presents a challenge for any model and suggests a re-examination of mechanisms which 
may cause very rapid variability would be worthwhile. In this paper I concentrate on 
how the observed variability time is related to the geometry and motion of the emis- 
sion region, and thus to the photon energy density in the emission region. The blazar 
emission mechanisms to be discussed include: a shock excited emission region, bent 
jets, a shock propagating along a jet containing a helical structure and illuminating 
parts of the helix by enhanced interactions/emission of radiation such that the emis- 
sion regions move along helical paths, and highly oblique conical shocks in the jet. 
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Together with geometry-specific time delays and variable Doppler boosting associated 
with relativistic motion of the emission region along a curved trajectory, it may well 
be possible to explain the observed flaring activity and high brightness temperatures. 
Another possibility briefly discussed in the context of bent jets and conical shocks is 
that jets may be fueled on an irregular time scale. 

The observed variability time At obs , and some assumed or estimated Doppler factor 
D is often used to estimate the jet-frame source radius, R' p=* DcAt ohs . The jet-frame 
photon energy density is then usually assumed to be U^ hot p=* L'/4ttR' 2 c, where V is 
the jet-frame bolometric luminosity given by V = D~ 4 4ird 2 L F, F being the observed 
bolometric flux, and being the luminosity distance. However, this approach can 
lead to large systematic errors in the jet-frame photon energy density. This is impor- 
tant because the energy density of photons of the low energy part of the SED may 
determine the energy losses of electrons, and the rate of up-scattering to 7-ray ener- 
gies in SSC models, and the rate of proton-photon collisions in hadronic models (see, 
e.g., Mannheim and Biermann 1989, Protheroe 1997, Mannheim et al. 2001, Miicke 
and Protheroe 2001, Miicke et al. 2002). If the emission region is optically thin, it 
may also have consequences for IDV sources as I will show that it is quite possible 
that the photon energy density responsible for the so-called Compton catastrophe may 
actually be lower than usually estimated. In the following sections I shall discuss how 
the variability time is related to the emission region geometry and intrinsic jet-frame 
variability, show how this can lead to large systematic errors in the jet-frame photon 
energy density, and discuss other scenarios which can lead to rapid variability and high 
fluxes. 

2 Relationship between variability time and emis- 
sion region geometry 

I shall discuss first a relativistic jet pointing at angle 9 with respect to the line of sight 
to the observer, in which the emission region is a cylinder of radius r and jet-frame 
length £' moving along the jet with the jet's Lorentz factor T = — (3 2 . To work 

out how the observed lightcurve depends on the emission region geometry and the 
duration of the emission in the jet frame (primed coordinates) we first consider in the 
observer frame (unprimed coordinates) the events corresponding to: (i) the emission 
from centre of cylinder at (t = = 0, x = x\ = 0, y = y x = 0, z = Z\ = 0) with this event 
defining the origin of coordinates in both frames, i.e. (£' = t[ — 0, x' = x[ — 0, y' = 
y[ — 0, z' = z[ = 0) ; (ii) the emission from some arbitrary point in the cylinder at a 
later time (t ii? Xa, ya, z^) ; (iii) the arrival at the telescope, located in the x-y plane, of 
the photon emitted in event (i); (iv) the arrival at the telescope of the photon emitted 
in event (ii) (see Fig. 1). 

In this paper, I shall assume that the distance between the AGN and the observer 
is very much larger than the dimensions of the emission region (R and £ in Fig. 1), 
and that we are not concerned with investigating variability on time scales very much 
shorter than R/c or £/c. Under this assumption, the time interval between arrival of 
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the two photons at the telescope is independent of distance to the AGN and is simply- 
given by 

^obs = (*iv - = (tii ~~ cos 9 sin 9^j . (1) 

The error made by using this approximation is ~ (R/ 'di)(R/ 'c) and is negligible com- 
pared with the variability time scale associated with the emission region geometry 
investigated in this paper which is ~ (?"/c). 

Lorentz transformation to the observer frame, gives t„ = r(t(; + /3xy), xy = r(x[ i + 
flt'x) and = y^. Noting that the Doppler factor is D = [T(l — PcosO)]^ 1 = T(l + 
(3 cos 9') , and using the Aberration formulae, 

r i 

costf = — (cos0' + /3), sin # = — sin 9' (2) 

cos 9' = r\D(cos0-/3), sm9' = Dsin9 } (3) 

Eq. 1 becomes 

Us = r^ + ^-r^ + ^^-Ccose' + ^-jZ—sine' (4) 

Dc Dc 

Dt ohs = 4- ^cos0'- ^sin0'. (5) 
c c 

Let us suppose that the cylindrical emission region emits radiation simultaneously 
and uniformly throughout its volume between times t' = —At' and t' = At', as mea- 
sured in the jet frame. While a simultaneous burst violates causality, it is nevertheless 
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a useful case to consider because it enables us to determine clearly the contribution 
of emission region geometry to the observed variability time. For 9' < n/2 the first 
photon to arrive would have been emitted at (t' = —At', x' = £'/2, y' — r, z' — 0) 
giving 

t^f = -At' - (£'/2c) cos 9' - (r/c) sin 9'. (6) 

Similarly the last photon to arrive would have been emitted at (f = At', x' = 
—£'/2, y' = — r, z' = 0) giving 

t^ = At' + (£'/2c) cos 9' + (r/c) sin 9'. (7) 

For 9' > 7r/2 the first photon to arrive would have been emitted at (t' = —At', x' = 
—£'/2, y' = r, z' = 0), and the last photon to arrive would have been emitted at 
(t' = At', x' = £'/2, y' = —r, z' = 0). Hence, if we define the observer frame duration 
of the burst as 2At obs then 

DAt ohs = At' + (£'/2c) | cos 0'\ + (r/c) sin 6' (8) 

and this is valid for all 8'. Note that Eq. 8 gives the usual formula, At G b s = At' /D, 
if the emission region is point-like (i.e. £' — r — 0). If one term in Eq. 8 dominates, 
DAt ohs gives one of: At', (£'/2c) \ cos6>'|, or (r/c) sin^'. 

I shall next consider the case of a cylindrical emission region in the jet being rapidly 
energized by a plane shock with jet-frame speed l3' shock c travelling along the jet, such 
that photons are emitted immediately after shock passage from a thin disk-like region 
immediately downstream of the shock. In this case, the location of the emitting disk 
is defined by x' = Pi^oc^ct' , and so the arrival times of photons at the telescope may be 
obtained from Eq. 5 

cDt ohs = x'(—^ costf') -y'sinO'. (9) 

V Pshock / 

For (l//9ghock ~ cos#') < 1 the first and last photons to be received would have been 
emitted at (x' = £'/2,y' — r,z' — 0) and (x' = —£'/2,y' = —r,z' = 0), respectively. 
However, for (l/f3' shock — cos 9') > 1 the first and last photons to be received would 
have been emitted at (x' = —£'/2,y' — r,z' — 0) and (x' = £'/2,y' = —r,z' = 0), 
respectively. Hence, the observer frame duration of the burst is 



DAt obs = (f /2c) 



-cos#' 



8' 

/^shock 

If 9' is small (9 is very small), one finds 

1 - B' 

A^shock 



+ (r/c) sin 9'. (10) 



At 



obs 



2Dc 



shock 



r9' 



and for a reasonable shock speed, e.g. f3' shock ~ ±0.5, 



£' r9' 

At* ~ PT1)^ + S (12) 
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where (3' shock ~ — 0.5 corresponds to a reverse shock. 

What we have learned from this discussion is that multiplying cAt Q bs by the Doppler 
factor might give the jet-frame intrinsic variability time or one of the dimensions of the 
emission region (possibly multiplied by some unknown factor). One is tempted to ask 
if R' m DcAt b s is the right dimension to put in Up hot = L' /4ttR' 2 c in order to estimate 
the jet-frame photon energy density. We shall discuss this point further in Section 4. 

3 Monte Carlo investigation of cAt b s , size, and D 

The Monte Carlo method allows the accurate calculation of expected light curves for 
any emission region geometry and intrinsic source variability. The emission region is 
modeled in the jet frame, and is represented by iV "particles", each of which emits pre- 
cisely one "photon" . The number density of the particles models the geometry of the 
emission region, and each particle emits its photon at a time determined by the emis- 
sion region geometry and variability model. The jet-frame 4-position of each photon 
emission event (^, x[, y[, z^), i = 1, 2, . . . N, is determined by model. The emission re- 
gion moves in the x-direction with Lorentz factor T, and the observer-frame 4-position 
of each photon emission event (ti, Xi, yi, Zi), % — 1, 2, . . . N, is obtained by Lorentz trans- 
formation. The arrival time of each photon is calculated for a given viewing angle 9 
using Eq. 1, and is binned to give the lightcurve. 

Two emission region geometries are considered: a solid cylinder, and a spheroidal 
3D gaussian. We also consider four intrinsic jet-frame time distributions: a simulta- 
neous burst (violates causality), a simultaneous gaussian pulse (violates causality) and 
simultaneous emission at plane shock. I shall investigate effects of varying the Doppler 
factor D, viewing angle 9, jet Lorentz factor T, and shock speed /9 s hock- Finally, the 
effect of smoothing due to acceleration/radiation time delays is discussed. 

3.1 Solid cylinder emission region geometry 

If the emission region is "solid" , i.e. the emissivity is constant inside the emission region 
and zero outside, the lightcurve may be peaked and have finite duration to reflect the 
sharp edges of the emission region, and the shape of the lightcurve will change with 
with viewing angle. The Monte Carlo results for the cylinder, when plotted such that 
the time of observation is divided by the expected variability time At obs given by Eqs. 8 
and 10, is shown in Fig. 2(a) for a simultaneous burst (At' = 0) , and in Fig. 2(b) for 
shock excitation. As can be seen, Eqs. 8 and 10 are verified by the Monte Carlo results. 
In both cases the shape of the lightcurve depends strongly on the viewing angle. 

3.2 Gaussian emission region geometries 

The density of "particles" representing the emission region is described by 




(13) 



6 



(a) Cylinder (M' = 0) 
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(b) Cylinder (/S' shock = + 0.5) 



DcAf ob ,=rsinfl' + 0.5Z'|l/|S , lllook -cosfl' 
r/(2r) = 0.1 

r=io 



r — fr 



•■'a 





obs 



Figure 2: Lightcurve from cylindrical emission region with ijlr = 0.1 moving with 
Lorentz factor T = 10 as observed at various viewing angles: (a) Simultaneous burst for 
& = 0° (solid), 36° (dotted), 72° (short-dashed), 108° (short-dashed, identical to 72°), 
144° (dotted, identical to 36°), & = 360° (long-dash, identical to 0°); (b) Excitation 
by a plane shock moving at speed /5g hock c = 0.5c for 9' = 0° (solid), 36° (dotted), 72° 
(short-dash), 108° (dot-dash), 144° (dot-dot-dot-dash), 9' = 360° (long-dash, identical 
to 0°). 

with a x = ay corresponding to a spherical gaussian distribution, a x < a r to an oblate 
spheroidal gaussian distribution, and a x > a r to an prolate spheroidal gaussian dis- 
tribution. For the case of a simultaneous gaussian pulse, the probability of emission 
at jet-frame time t' to (t' + dt') is p(t')dt', where pit') = exp(— t' 2 /2a t 2 ) / V^cft is 
independent of position. 

If the emission region has a gaussian shape the lightcurve will be smooth, and in 
many cases will also have a gaussian shape. This is true, for example, for the cases 
of a simultaneous gaussian pulse and excitation by a plane shock. The width of the 
lightcurve will depend on viewing angle in approximately the same way as for equivalent 
cylindrical emission volume. For example, r and £'/2 of the solid cylinder should be 
related to ay and a x , respectively, for the case of the spheroidal gaussian density. For 
the case of a simultaneous gaussian burst I find the standard deviation to be given by 

a = [(a r sin 9' / cf + (a x cos 9'/cf + a^/D. (14) 

Note the similarity to Eq. 8 except that the terms are added in quadrature as the 
standard deviation is required instead of the maximum duration of the pulse which is, 
theoretically, infinite. 

In the case of excitation by a plane shock wave with speed (3 s hockC, the emission 
is from the plane where the shock cuts the spherical gaussian density. The emission 
region then has a surface density of emitting particles which is a two-dimensional 
gaussian, and the lightcurve reflects this distribution, and so is also gaussian. The 
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standard deviation of the gaussian lightcurve depends on the viewing angle as a result 
of projection effects, and I find the standard deviation to be 



a 



1 

cD 



(a r sm6') 2 + (a x y 



cos 9' 



shock 



1/2 



(15) 



Lightcurves for T = 10, f3' shock = 0.5 and various viewing angles are plotted in Fig. 3 
and seen to lie on top of each other when plotted in units of a, and to be a normal 
distribution. Equation 15 is valid for both forward and reverse shocks. For reverse 
shocks, /Sshock would be negative and give rise to a broader lightcurves (larger a). 



spherical gaussian density 




Figure 3: Lightcurve due to excitation by a plane shock of a spherical density moving 
with Lorentz factor T = 10 for f3' shock = 0.5 and various viewing angles shown is seen 
to be a normal distribution. 



3.3 Acceleration/radiation delays 

Instantaneous excitation by a shock wave is not a realistic approximation unless the 
time-scales for particle acceleration and radiation (cooling) are very short compared to 
the transit time of the shock through the emission region. If this is not the case, the 
lightcurve would be broadened and smoothed to reflect the time delays associated with 
particle acceleration and radiation. An equation of the form dN'/dt' = Q'(t') —N'/t[ oss 
might describe the time-evolution of the number of particles, N', radiating photons 
at the frequency corresponding to that observed, with Q'(t') being the source term, 
and t( oss representing the time-scale for particle losses (or time-scale over which the 
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radiation is emitted). Then a simple smoothing function of the form 
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can be used for this purpose, where t' dei represents the delay between shock passage and 
emission by the radiating particles, At gain represents the the duration of the acceleration 
following shock passage. In the case of a pre-existing population of thermal electrons 
shock heating is essentially instantaneous (At gain = 0), and t( oss represents the time- 
scale for cooling by thermal bremsstrahlung radiation. Similarly, for a pre-existing 
population of relativistic electrons, the passage of the shock will essentially instantly 
increase the magnetic field (At gain = 0) and t( oss would represent the time-scale for 
energy losses by synchrotron emission. If particle acceleration is required then At' 



gam 



> 



0. 

Fig. 4 shows the effect of acceleration/radiation delays on the lightcurve due to 
a plane shock (fl' shock = 0-3) exciting an oblate gaussian spheroidal emission region 
(o x jo r = 0.2) with acceleration/radiation delays described by At gain = t| oss for t| oss = 
O.lay/c, ay/c and 10oy/c, and viewed at angle 9' = 5°. In the case of no accelera- 
tion/radiation delays, the lightcurves would simply be normal distributions centred on 
^obs = 0. Results are shown for the three t( oss - values, and t G bs is plotted in terms of cr to t 
defined by 
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(17) 



where a is given by Eqn. 15, such that a tot gives a crude measure of the expected 
duration of the lightcurve. These distributions can be obtained simply by convolution 
of a normal distribution with de i (Eq. 16), taking account of the fact that jet-frame 
times enter in Eq. 16, that for any point co-moving with the jet At obs = At'/D, and 
that t bs is plotted in units of a tot , 



( tobs \ 
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0dei(iobs-D - 1) exp 



2a 2 



dt. 



(18) 



The solid curves, which agree with the Monte Carlo results (histograms), are obtained 
from Eqn. 18. As we see, if the time-scales for the acceleration/radiation process (At gain 
and t( oss ) are much less than the time-scale a' associated with the shock passage and 
dimensions of the emission region then the light-curve will be symmetrical, and in the 
case of a gausian spheroidal emission region will be a gaussian distribution with stan- 
dard deviation a (leftmost histogram). If the time-scales for the acceleration/radiation 
process are much larger than a', then the light-curve will reflect that of the accel- 
eration/radiation process, i.e. Eqn. 16 (rightmost histogram). The middle histogram 
shows an intermediate case. 

Lightcurves of AGN show many different features and forms of variability. One 
example which appears to show variations reflecting the acceleration/radiation process 
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Figure 4: Excitation of an oblate spheroidal gaussian density (a x j o r = 0.2) moving with 
Lorentz factor T = 10 by a plane shock with jet-frame speed /9ghock c = 0.3c, followed 
by particle acceleration and radiation represented by Eq. 16 with Atg ain = t( oss for 
t[ oss = O.lov/c, cr r /c and 10av/c (as indicated) for viewing angle 9 = 5°. Solid curves 
give the result of convolution (Eqn. 18). 



is 3C 454.3. In Fig. 5, I show the lightcurve of 3C 454.3 at 37 GHz obtained with the 
Metsahovi and Crimea telescopes over ten years (Salonen et al. 1987, Terasranta et 
al. 1992). The flares appear non-symmetrical and have a shape similar to the right- 
most histogram in Fig. 4 (Atg ain = t( oss 3> a). As an example, I have constructed 
a reasonably well-fitting lightcurve (solid curve) from a number of flares of the form 
given by Eqn. 16 with At gain = ti oss = 0.75 yr (dotted curves) plus a background flux 
density of 3 Jy. One interpretation of these data would then be that the radio emission 
region in 3C 454.3, being modelled by the solid curve, has dimensions much less than 
c(0.75 yr)D = 0.1 QD pc, and that the energy-loss time-scale of the radiating electrons 
is ~(0.75 yr)L>. 

Another example is the X-ray lightcurve of 1ES1959+65 from ARGOS/USA and 
RXTE/ASM (Giebels et al. 2002) shown in Fig. 6. In this case, the individual flares 
appear to be symmetrical in time, and have a roughly gaussian shape. I have con- 
structed a reasonably well-fitting lightcurve (solid curve) from a number of flares of 
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Figure 5: Lightcurve of 3C 454.3 at 37 GHz (Salonen et al. 1987, Terasranta et al. 1992). 
The fit (solid curve) comprises a background level (horizontal dotted line) and several 
flare components (dotted curves) with individual flare lightcurves given by Eq. 16 with 
At gain = ti oss = 0.75 yr. 



gaussian form with standard deviations a = 2.5 d (dotted curves) plus a background 
flux density of 2 mCrab. One could interpret this as indicating that the energy-loss 
time-scale of the radiating electrons is much less than 2.5.D d, and that the emission 
region dimensions are ~ c2.5D d= 2 x 10~ 3 D pc. 

4 Dependence of energy density on dimensions of 
the blob 

It is important to know the dimensions of the emission region for several reasons: (i) in 
some hadronic models the synchrotron photons are targets for photoproduction, (ii) in 
all models the synchrotron photons are targets for photon-photon pair production by 
7-rays, (iii) in SSC models the synchrotron photon energy density determines Compton 
scattering, and (iv) knowing the systematic errors on photon energy density may help 
understand the so-called "Compton catastrophe" in IDV sources which have apparent 
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Figure 6: X-ray lightcurve of 1ES 1959+65 (Giebels et al. 2002). The fit (solid curve) 
comprises a background level (horizontal dotted line) and several flare components 
(dotted curves) with individual gaussian lightcurves with standard deviations a = 2.5 d. 



brightness temperatures well in excess of the limit T B < 10 K imposed by Compton 
scattering (Kellermann and Pauliny-Toth 1969, Kardashev 2000) when the photon 
energy density in the emission region reaches the energy density in the magnetic field. 
One extreme example is PKS 0405-385 (Kedziora-Chudczer et al. 1997) which has 
T B > 5 x 10 14 K after correcting for interstellar scintillation (see also Walker 1998) 
requiring a Doppler factor of D = 10 3 to satisfy the brightness temperature limit. 

To illustrate how critically the energy density depends on the geometry, I shall 
consider the case of the jet-frame emissivity following a spheroidal gaussian density. 
Provided the emission is optically thin, as is almost certainly true for the optical-X-ray 
synchrotron hump in the SED of blazars, then it is straightforward to calculate the 
average energy density from the emission region geometry and the luminosity. 

Assuming that the jet-frame luminosity, L', is constant, we can estimate the average 
jet-frame photon energy density given the Doppler factor and emission region geometry. 
The simplest way of doing this, for any emission region geometry, is to use the Monte 
Carlo method to place N points at positions r iy i = 1, . . . N, distributed according to 
the emission region geometry, and to give each point a luminosity L'/N. Then at point 
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fj the energy density is 

^)^I>- (19) 
where = \ri — fj\. Averaging over the emission region distribution we obtain 

Tl N 

Khot>~^EE^ 2 - (20) 

i=i j^t 

The result giving the jet-frame average photon energy density versus a x is shown as 
the thick solid curve in Fig. 7 for the case of a gaussian spheroidal density. 
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Figure 7: Jet-frame average photon energy density versus o x . Thick solid curve shows 
actual density calculated using the Monte Carlo method. "Observed" density is shown 
by the shaded areas for a t = (top shaded area), a r , 4<r r , 8<r r , and 16cr r (bottom 
shaded area). See text for further details. 



In Fig. 7 we also plot the jet-frame average photon energy density that would be 
inferred if we assumed that the observed variability time scale a and an assumed or 
estimated Doppler factor D gave the jet-frame radius of a spherical emission region, i.e. 
R' = aDc. This "observed" jet-frame average photon energy density is simply given 
by 

u 

(^hot.obs) = 47r((JjDc) 2 c - ( 21 ) 

Note that a (given by Eq. 14) depends on 6', a t and o x so that the "observed" jet- 
frame average photon energy density depends also on 9', a t and a x as well as D. 
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The figure shows (^phot,obs) plotted against a x and gives the range due to variation 
in 9' (shaded) for various a t . The upper bound in each case gives the result for 9' = 
or 7r and the lower bound is for 9' = n/2. We see that the "observed" value, 
(^phot.obs) can be several orders of magnitude higher or lower than ({/ phot ) if the emission 
region is different from a sphere, or if the intrinsic variability time a t is not small. 
For example, take the case of a t <C cr r /c, if a x = 10~ 2 o" r and 9' = (9 = 0) then 
(^hot,obs) ~ 10 3 (^ hot ), whereas if a x = 10V r and & = 90° (9 = 5.74° for T = 10) 
then (C^phot,obs) ~ 10 2 (^ P hot)- The above result has clear implications for both leptonic 
and hadronic models of AGN in which photons of the low-energy peak of the SED 
provide target photons for inverse- Compton scattering by electrons (leptonic models) 
or pion photoproduction by protons (hadronic models). Using the observed variability 
time together with assumed or estimated Doppler factor to estimate the emission region 
radius R' can clearly lead to large errors in (£/ phot ). 

Although the radio emission in IDV sources is usually assumed to be optically thick, 
if this is not the case then the above result may also have implications for IDV sources 
as the photon energy density responsible for causing the brightness temperature limit 
may actually be a few orders of magnitude lower than estimated on the basis of time- 
variability, and in that case much lower Doppler factors would be required to avoid the 
Compton catastrophe. I explore this further in a separate paper (Protheroe 2002). 

5 Variability due to photon pile-up in observation 
time 

The simplest example of photon pile-up in observation time is a bent jet. Jets may 
bend if they pass through a stratified cold and high density region (Mendoza and 
Longair 2001). We approximate the trajectory of an emission region moving with 
speed (3c = (1 — l/T 2 ) -1 / 2 along a bent section of jet by motion around a section of 
a circle of radius r in the x-y plane: (t, x = rcos(ut), y = rsin(cjt), z = 0), where 
uj = (3c/ r. For an observer in the x-z plane at angle r\ to the rr-axis, the observation 
time is then given by 

tobs = t — (r/c) cos(cjt) cost] (22) 

and the Doppler factor is 

D = {T {I + /3 cos rising)]}- 1 . (23) 

The Doppler factor raised to the 4th power is plotted against observation time for 
T = 10 and various observation angles in Fig. 8. I have plotted D 4 as it is appropriate 
for bolometric flux from a moving isotropic source; it also applies to the specific flux 
F u for F v oc v~ x . We see that even for modest observation angles the lightcurve 
is strongly peaked, essentially a delta function when the emission region direction is 
closest to the line of sight. Of course the finite size of any emission region will broaden 
the distribution. For example, for a lab-frame emission region length along the jet £, 
the burst would have duration At obs ~ £/c. 
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Figure 8: Doppler factor raised to 4th power representing bolometric lightcurve due to 
emission region moving along a bent jet. The jet Lorentz factor is T = 10, the bent 
jet is approximated by part of a circle around which jet plasma moves with angular 
velocity u>, and the emission region is observed at angle r\ < 0.3° (solid curve), rj = 27° 
(chain curve), r] = 81° (dashed curve). 



5.1 Helical jet structures 

VLBI observations show that helical jets or helical structures in jets may be fairly com- 
mon in AGN (Rantakyro et al. 1998), and theoretical studies have shown that wave-like 
helical structures can occur as a result of jet precession (Hardee 2000). Several papers 
discuss helical jet models or the application of helical models to specific sources (Ca- 
menzind 1986, Rosen 1990, Tateyama et al. 1998, Qian et al. 1992, Schramm et al. 
1993, Steffen et al. 1995, Villata and Raiteri 1999). Certainly helical jets or struc- 
tures would be important in determining the lightcurve of 7-ray and neutrino emission 
from blazars, and various suggestions have been made about the mechanisms involved 
(Despringre and Fraix-Burnet 1997, Marcowith et al. 1995). 

I consider a filamentary helical structure embedded in the jet with Lorentz factor T 
whose axis coincides with the jet axis and is excited by a plane shock travelling along 
the jet with jet-frame speed f3' shod!i c. The helical structure could be, for example, a 
flux tube containing a relatively high magnetic field, or a tube of high plasma density 
arising from a density perturbation in the plasma entering the jet. Helical magnetic 
fields with an Archemedian spiral topology similar to the "Parker spiral" field of the 
heliosphere may well be expected in AGN jets. 

The excitation of a filamentary helical structure is described by the jet-frame 4- 
vector (t' , x' , y' , z') = (X(fi/ /3' shock c, A0, rcos0, rsin0), where r is the radius of the 
cylinder containing the helix, A is the helix wavelength, A" = A/27T, = f3' shock ct' / X , 
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and the jet is pointing in the x-direction. Lorentz transformation to the galaxy-frame 
gives (t,x,y,z) = [rt'(l+PP' shock ) : Tct'((3+(3' shock ), r cos(f3' shock ct' / X) , rsm((3' shock ct'/X)}. 
The galaxy-frame speed of the location of the excited part of the helix can exceed c, 
but this does not violate causality as no particles or information propagates at this 
pattern speed which is 



/^shock + P 
^pattern 1 R , R 
1 "T PshockP 



/ rfl' 

A^shock 



*r(/4 ock + /3). 



(24) 



Observation in the x-z plane at angle 9 to the jet axis (rr-axis) gives 

t obs = r^i + /?/4 ock - (/4 oc ^ (25) 

^ = T[l + PP' shock - cos9(/3' shock + (3)} +rsm9(P' shock /X) cos(/4 ock ct'/X).(26) 

Whenever dt \, s /dt' = the lightcurve will have a cusp. If cusps are possible, they will 
occur at times corresponding to 

, -T[l+(3(3' shock -cos9(P' shock + (3)] 
" rsin0(/4 ock A) (27) 

provided the model parameters give cos0 in the range — 1 < cos0 < 1 (cos0 depends 
on the helix geometry, shock speed, jet Lorentz factor and viewing angle). If cos = ±1 
one cusp per helix wavelength will occur, and if — 1 < cos0 < 1 multiple cusps occur, 
otherwise no cusps are present in the light curve. However, if cos0 is close to ±1 the 
lightcurve will be peaked. This is illustrated in Fig. 9 for T = 10, /3g hock = 0.5, and 
A = r. Fig. 9(a) shows t' plotted against t Q b s for five viewing angles, and Fig. 9(b)-(f) 
shows the resulting lightcurve for each of the five viewing angles (the corresponding 
Doppler factors are also given). 

Apart from the periodicity, the lightcurves shown in Fig. 9 are reminiscent of those 
of blazars. All of these lightcurves correspond to instantaneous emission from the point 
on the helical filament at the time of excitation. The cusps and peaks would in reality 
be smoothed to some extent by the finite width of any helical structure as well as 
by any delays associated with acceleration and radiation time scales. Furthermore, 
as the shock weakens while propagating down the jet the successive peaks/cusps in 
the lightcurve would decrease in height. One possible scenario could be that there 
is a succession of shocks at random intervals propagate down the jet, and because 
individual shocks weaken as they propagate, each shock would cause only one or two 
peaks (due to only one or two cycles of the helix). In this case, it may not be possible 
to distinguish between a simple bent structure in the jet and a helical structure - if 
appropriately aligned, a simple bent structure would cause a cusp in the lightcurve in 
exactly the same way as a helical structure. An alternative to a helical filamentary 
structure within the jet is a helical jet which would produce a qualitatively-similar 
lightcurve, but would have the Doppler factor varying with position along the jet as 
the viewing angle relative to the local jet direction changes. 
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5.2 Conical shocks 



Lind and Bladford (1985) have considered the possibility that hotspots seen in VLBI 
images of radio jets may actually be relativistically moving conical shocks. They ap- 
plied the relativistic shock jump conditions to a plane parallel flow entering a forward 
conical shock with cone angle 77 to find the angle ( < rj with at which the flow initially 
diverges with respect to the jet axis. Defining the downstream region between the 
cones with angles ( and 77 to be the emission region, and taking account of Doppler 
boosting, they model the brightness distribution to simulate VLBI images. 

The generation of conical shock structures is often seen to occur in simulations of 
relativistic jets after the introduction of a fast perturbation (Bowman 1994, Gomez 
et al. 1997). These structures, which are typically an alternating sequence of forward 
and reverse conical shocks, are stationary or slowly-moving in the galaxy frame, and 
have cone angles 77 ~ 1/T. Salvati et al. (1998) discuss emission by a conical shock. 
They consider the case of a density perturbation, confined to thin flat disk, travelling 
relativistically along the jet causing particle acceleration and emission where the disk 
cuts a stationary forward conical shock such as that illustrated in Fig. 10. 

Salvati et al. (1998) assume that the time-scales for acceleration and emission are 
negligible. They show that their model can lead to highly peaked light curves whose 
shape depends on the viewing angle, and used these light curves to fit Markarian 421 
TeV flare data. For the same input, and assuming there is no Doppler boosting of 
the emitted radiation, I am able to reproduce exactly their figure 2 which shows the 
observed flux for various viewing angles. However, bearing in mind that the shock is 
stationary and that the pile-up in observing time is already included, the bolometric 
flux emitted by part of the downstream flow will be Doppler boosted by -Df ocal where 
Aocai = [r<f(l — Aicos^iocai)] -1 where T d is the Lorentz factor of the downstream 
flow, and 6>i ocal is the viewing angle with respect to the line of sight and the local 
downstream flow direction which varies around the shock as indicated in Fig. 10 (0i OC ai 
ranges between 9 a and 6b). To obtain the Lorentz factor of the downstream flow, it 
is easiest to Lorentz transform in a direction parallel to the shock plane to a frame 
in which the flow is normal to the shock. For the case of cone angle 77 = sin -1 (1/T) 
assumed by Salvati et al. (1998), i.e. an oblique shock at angle 77 to the upstream flow, 
and using the relativistic equation of state, I find that for T = 10, 77 ~ 7.25° and 
( 1.89° and that the Lorentz factor of the downstream flow is related to that of 
the upstream flow by T d « 0.801784r u . My result for the observed flux for the same 
input as Salvati et al., but including the Doppler boosting taking into account the 
local downstream flow directions, is given in Fig. 11(a) and shows that the inclusion 
of Doppler boosting causes the peak at 9 = to be higher than that at — 0.977, 
the opposite to that found by Salvati et al. However, since Salvati et al. (1998) used 
one viewing angle in their fits to Markarian 421 TeV flare data (their figure 3), and 
because the divergence of the downstream flow is rather small, their fits are still valid 
and their model remains an interesting mechanism for flare production. The same 
authors (Spada et al. 1999) have applied their model to IDV sources and are able to 
explain brightness temperatures up to 3 x 10 17 K. 
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I wish to extend the work of Salvati et al. (1998) by including the reverse shocks, 
and ultimately a sequence of stationary reverse and forward shocks. In Fig. 11(b) 
I show the lightcurve for a reverse conical shock having identical parameters as the 
forward shock already discussed. Because the divergence of the downstream flow is 
rather small, the lightcurve of the reverse shock is approximately just the lightcurve of 
the forward shock reflected about t G b s = 0. 

I shall consider next the case of a sequence of stationary reverse and forward conical 
shocks. Although in real AGN jets it may be possible that the downstream flow is re- 
accelerated to near the original value (depending on conditions external to the jet, 
and the jet production mechanism), I shall assume that each successive conical shock 
causes a reduction in the jet Lorentz factor by Td/T u = 0.801784. I shall make an 
additional approximation that the divergence/convergence of the flow caused by the 
conical shocks can be neglected, and that the jet flow is always parallel to the jet axis. 
The lightcurve due to a thin, initially flat, disk travelling relativistically along the jet 
is then calculated by the Monte Carlo method, in which "particles" are placed over the 
surfaces of all of the cones. Each cone has the same number of particles, and they would 
appear uniformly distributed across the cross section of the jet when viewed along the 
jet axis. A side view showing the location of these particles is given in Fig. 12. An 
initially flat thin disk is launched along the jet with initial Lorentz factor T = 10, and 
each time a part of the disk crosses one of the shocks the Lorentz factor of that part 
of the disk drops appropriately, distorting the disk. The time at which each "particle" 
emits its photon is determined by the time at which the (distorted) disk reaches the 
"particle". The resulting lightcurve for various viewing angles is shown in Fig. 13(a) 
where the emission is boosted using the Doppler factor corresponding to the Lorentz 
factor of the flow immediately downstream of each shock. If the Lorentz factor of 
the jet decreases at each shock crossing as assumed here then, because the Doppler 
factor is also reduced at each shock crossing, only the first few conical shocks would be 
prominent in the lightcurve resulting from a single jet perturbation. This is shown in 
Fig 13(b) which gives the lightcurve from Fig 13(a) for 9 = 9° on a linear scale. 

As can be seen, the lightcurve is quasi-periodic. It would be strictly periodic if the 
downstream flow velocity were identical to the upstream flow velocity. This might be 
the case if the re-collimmation of the diverging flow from the conical shocks results in 
restoration of the flow Lorentz factor to near the upstream value, and in this case the 
flux from successive cycles would be at about the same level, as observed in Markarian 
501 flares. The two flow velocities would also be roughly the same if the initial jet 
Lorentz factor were much higher than T = 10 used in Fig. 13(a) such that the intervals 
between flares due to a pair of conical shocks was roughly the same. However, in 
this case the flux of successive flares would diminish as the Doppler factor decreases. 
Turning this argument around, we may have a method of determining the minimum 
jet Lorentz factor from observational data, e.g. the ~ 23 day periodicity in the 1997 
Markarian 501 data (Protheroe et al. 1998, Hayashida et al. 1998) may be used to put 
a lower limit on T. 

Fig. 14 shows an example lightcurve due to multiple jet perturbations. In this case, 
I have used an exponential distribution of times between the injection into the jet of 
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a density perturbation with a mean interval of r/c. I have also used an exponential 
distribution for the strength of each perturbation. The resulting lightcurve looks, at 
least qualitatively, as good as any model for the variability of fluxes from blazars. Note 
that even though there is no strict periodicity resulting from the excitation of a sequence 
of shocks in this case, each flare episode has two strong peaks due to (for this viewing 
angle) the apex of the first reverse shock cone and the apex of the first forward shock 
cone, and in each flare the separation between these two peaks is identical. Fourier 
analysis would therefore show a strong peak at the frequency corresponding to the time 
interval between these two peaks in a single flare. 

6 Conclusion 

Many factors can influence observed variability time. The connection between At Q b s , 
Doppler factor and emission region geometry is non-trivial, and so measuring At Q b s 
may give, at best, one of the dimensions of the emission region. Using R' = cAt Q b s D 
and Up hot = L'/4ttR' 2 c may then lead to over-estimation or under-estimtion of the 
jet-frame photon energy density by orders of magnitude. This is clearly of importance 
in any AGN model in which the low energy photons produced in the jet are targets for 
interaction of high energy particles or radiation, such as in SSC models and hadronic 
blazar models. Although not discussed in detail in the present paper, the escape of 
7-rays from the emission region depends on the optical depth to photon-photon pair 
production interactions. This optical depth can be uncertain by orders of magnitude 
in the same way as the photon energy density, and will also depend on viewing angle. 
One must therefore be careful when using the observation of apparently unattenuated 
gamma-rays, and an observed variability time, to place limits on the Doppler factor. 

The uncertainty in the jet-frame photon energy density discussed in this paper 
may also have implications for the high brightness temperature/Compton catastrophe 
problem of IDV sources. In this case, it is the energy density of target photons which 
limits the brightness temperature through the competition of inverse-Compton scat- 
tering with synchrotron radiation, and the target photon energy density may actually 
be lower than estimated if the emission region is non-spherical. 

If the jet is bent or helical, or has some other favoured geometry (e.g. conical 
shocks) cusps in t' vs. t bs, and/or a varying Doppler factor may cause narrow peaks in 
the observed lightcurve irrespective of other factors. Distinguishing between these cases 
from the observed lightcurve alone is likely to be difficult. One way of distinguishing 
whether a flare is due to (i) an emission region moving around a bent or helical path, 
or (ii) a shock is exciting a curved, conical or helical structure within a jet, is that the 
in the first case the flare is caused by a change in viewing angle with respect to the 
motion of the emission region leading to a change in Doppler factor, whereas in the 
second case the flare is due to a pile-up in observation times with no change in Doppler 
factor. Hence, in case (i) not only will the observed flux increase during a flare, but 
the photon energies also increase - increase in (^i^) pea k by factor x A accompanied by 
shift in fpeak by factor x (x is ratio of final to initial Doppler factor). In case (ii), 
however, since there is no change in Doppler factor there should be no shift in u pcak 
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accompanying an increase in (z/Fj y ) pea k- 

Distinguishing between the excitation of a conical and a helical structure by a plane 
shock would be almost impossible - note the qualitative similarities between lightcurves 
depicted in Figs. 9(d) and 13(b). Relativistic jet simulations (Bowman 1994, Gomez et 
al. 1997) do show the presence of conical shocks, and these shocks appear after a large 
perturbation (e.g. from 4 to 10 in the simulations of Gomez et al. 1997) in Lorentz 
factor of the matter entering the jet through the nozzle. This can result in a sequence 
of quasi-stationary to superluminal reverse and forward conical shocks extending from 
the nozzle to the perturbation as it moves along the jet (Agudo et al. 2001). For 
the conical shock model (Salvati et al. 1998) discussed here to work, a subsequent 
perturbation would need to result in a plane shock, or plane thin perturbation of some 
kind, which could travel along the jet and excite the pre-existing conical shocks. As far 
as I am aware, whether or not this could occur has not been demonstrated. A similar 
uncertainty hangs over whether or not helical jet structures, which may themselves 
be shocks with a twisted ribbon topology, resulting from a perturbation entering the 
jet through the nozzle, could subsequently be excited by the passage of a plane shock 
or perturbation. Nevertheless, in both cases, if the viewing angle is favourable pile- 
ups in t bs, and hence flares, could occur simply as a result of the motion of the 
conical or helical patterns which may themselves be sites of enhanced emission. Note 
that in recent 3D relativistic jet simulations, the introduction of a 1 percent helical 
velocity perturbation at the nozzle results in a helical pattern propagating along the 
jet at nearly the beam speed (Aloy et al. 1999), and that the conical shocks resulting 
from a perturbation in jet Lorentz factor can range from being quasi- stationary to 
superluminal (Agudo et al. 2001). 

In conclusion, in models for flaring in AGN in which the emission comes from a 
localized region (blob) co-moving with the jet, time variability is non-trivial to interpret 
in terms of emission region geometry and Doppler factor. A further complication is 
that flaring may arise instead due to curved or helical motion of a blob, even if the 
emission is constant in the istantaneous rest frame of the blob. In this case, apparent 
flaring is due to the change in viewing angle, and hence Doppler factor. Similarly, if 
the viewing angle is favourable, relativistic motion of curved or helical filaments or 
surfaces can lead to observation of flares. Excitation of curved or helical jet structures 
by shocks or perturbations can also lead to pile-ups in t b s , and hence large apparent 
increases in flux. Observations of time variability in AGN is therefore non-trivial to 
interpret and may lead to large systematic errors in estimated jet-frame photon energy 
density, Doppler factor and the physical parameters of the emission region. 
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Figure 9: Shock excited helical jet structure for T = 10, f3' shock = 0.5, A = r. (a) t' vs. 
t bs for 9 = 0.3° (solid curve), 1° (dotted curve), 3° (short dashed curve), 9° (dot-dash 
curve), and 27° (dot- dot- dot- dash curve). Resulting lightcurves for (b) 6 = 0.3°, (c) 
1°, (d) 3°, (e) 9°, and (f) 27°. 
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Figure 10: Side view of the jet showing the location of a stationary forward conical 
shock, and showing the downstream flow diverging from the jet axis at angle ( (not 
to scale). For a viewing angle 9 with respect to the jet axis, the angle between the 
downstream flow direction and the line of sight varies between 9 a and 9b- 




Figure 11: Lightcurve of (a) a stationary forward conical shock and (b) a stationary 
reverse conical shock locally excited by a thin density perturbation travelling along 
the jet with Lorentz factor T = 10 as viewed at angle 9 with respect to the jet axis. 
The cone angle is rj = sin _1 (l/r) and results are given for 9/rj = 0.001 (long dashed 
curves), 0.1 (dot-dot-dot-dash curves), 0.3 (dot-dash curves), 0.5 (short dashed curves), 
0.7 (dotted curves), 0.9 (solid curves). The Doppler factor is calculated using the local 
downstream flow velocity. 



24 




Figure 12: Side view of the jet showing the location of "particles" used in the Monte 
Carlo Method. The location of the conical shocks is clearly evident. 
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Figure 13: Lightcurve due to excitation of the sequence of stationary conical conical 
shocks discussed in the text for an initial jet Lorentz factor T = 10 (a) as seen at 
viewing angles 6 = 0° (solid curve), 3° (dotted curve), 6° (short dashes), 9° (dot-dash) 
and 12° (dot-dot-dot-dash) with respect to the jet axis, and (b) as seen at viewing 
angles 9 = 9° and plotted on a linear scale. 
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Figure 14: Lightcurve due to repeated excitation of the sequence of stationary conical 
conical shocks discussed in the text for an initial jet Lorentz factor T = 10 as seen 
at viewing angles 6 = 9°. The time interval between perturbations entering the jet 
is sampled from an exponential distribution with mean r/c, and the strength of the 
perturbation is also sampled from an exponential distribution. 
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